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A number of LHC resonance search channels display an excess in the invariant mass 
region of 1.8-2.0 TeV. Among them is a 3.4 a excess in the fully hadronic decay of 
a pair of Standard Model electroweak gauge bosons, in addition to potential signals 
in the HW and dijet final states. We perform a model-independent cross-section fit 
to the results of all ATLAS and CMS searches sensitive to these final states. We 
then interpret these results in the context of the Left-Right Symmetric Model, based 
on the extended gauge group SU{2)l x SU(2)r x U( 1)', and show that a heavy 
right-handed gauge boson Wr can naturally explain the current measurements with 
just a single coupling qr ~ 0.4. In addition, we discuss a possible connection to dark 
matter. 
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I. INTRODUCTION 

A number of searches for narrow resonances at the LHC now display peaks in the region 
of 1.8-2.0 TeV. Most prominently, a recent ATLAS search [I] for a resonance that decays 
to a pair of standard model (SM) gauge bosons contains a local excess of 3.4cr (2.5 a global) 
in the WZ final state at approximately 2 TeV. Since the search is fully hadronic, there is 
only a limited ability to accurately distinguish gauge bosons, and consequently many of the 
events can also be interpreted as a ZZ or WW resonance, leading to excesses of 2.9 a and 
2.6 a in these channels respectively. 

Interest is further piqued when one analyses other resonance searches in a similar mass 
range. At approximately 1.8 TeV, both CMS [2] and ATLAS [3] observe excesses in the dijet 
distributions with a significance of 2.2a and la, respectively. In addition, a CMS search 
for resonant HW production [Hj shows a 2.1 a excess, and another CMS search [5j| for a pair 
of vector bosons, but this time with a leptonically tagged Z, finds a 1.5 a excess, both at 
approximately 1.8 TeV. With the exception of the three ATLAS selections in gauge boson 
pair production, all these possible signals are completely independent. 

We analyze these excesses and perform a general cross-section fit to the data. It is of 
course important to not only look at possible signals in the distributions, and for that 
reason we include all relevant searches that may be sensitive to the same final states as 
those showing anomalies. These include diboson analyses from both ATLAS [EHE] and CMS 
omj. In addition, many models that explain an excess of events in the clijet distribution 
with a charged mediator will also lead to a peak in a tb resonance search. Consequently, 
we also include the most sensitive searches for this particular final state from both ATLAS 
[12 , 13] and CMS [H] in our study. 

Combining all searches, we provide best fit signal cross sections for each of the final states 
analyzed in order to guide the model building process. In particular we find that the vector 
boson pair production searches are best described by a, WZ or ZZ final state, while WW is 
disfavoured via limits from semi-leptonic searches. In the associated Higgs production fit, we 
find preference for the HW final state, since an excess here is only seen in the single-lepton 
analysis. Finally, there is good agreement for a dijet signal between ATLAS and CMS, but 
nothing has been observed so far in the tb final state. 

In order to explain these excesses we focus on the so-called Left-Right Symmetric Model 
(LRM) pT5LIT9] . The LRM is based on the low-energy gauge group SU(2) L x SU(2) R x U(l)' 
that can arise for example from an S'O(IO) or E 6 GUT [20, 2T]. Since there is a new SU(2) 
gauge group in addition to the SM, the spectrum now contains both a Z R and charged Wf t 
bosons. In general the model predicts mz R > m w ±i which is why we explain the various 

| ^ 

signals via resonant W^ production and decay. By extending our cross-section fit to the 
LRM parameters, we find a region that can explain the excesses while avoiding all current 
exclusion bounds. We then examine where our fit suggests the mass of the Z R should be 
and the discovery potential of the LHC for this state. 

In a next step we also analyze the question of how such a model may be able to simul¬ 
taneously explain dark matter (DM). If the W R is to mediate DM annihilation in the early 
universe the simplest scenario we can find is to include a new charged and neutral particle 
that would be relatively close in mass in order to co-annihilate effectively. Alternatively, we 
offer the concrete example that DM annihilation is mediated through the Zr and the dark 
matter candidate is the neutral component of a new fermionic doublet. In this case, we see 
that resonant annihilation is required for the correct relic density and thus we predict the 
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DM mass is rv_/ • 

This paper begins in Sec. |TT] with a review of the different experimental studies under 
consideration and generic cross-section fits. In Sec. Ill we interpret the results in the LRM, 
followed by a discussion of a potential link to dark matter in Sec. IV Finally, we give our 
conclusions in Sec. El 


II. EXPERIMENTAL MEASUREMENTS AND CROSS-SECTION FITS 


In this section we give an overview over the ATLAS and CMS searches in which resonances 
in the mass range 1.8-2.0 TeV have been observed, or which provide constraints in the same 
or closely related final states. It is important to note that due to resolution, jet energy scale 
uncertainties and the limited number of events currently present in the excesses, peaks at 
1.8 -2.0 TeV can easily be compatible. Even though the most significant excess is observed 
in the ATLAS diboson search [I] at an invariant mass of 2.0 TeV, the majority of excesses 
are reported closer to an invariant mass of approximately 1.8 TeV. This is why we assume 
a resonance in the vicinity of 1.8-1.9 TeV in this analysis. 

We analyze the results of these searches by performing cross-section fits in each channel 
individually. For this, we sum bins around this mass region according to the experimental 
width and jet energy scale of the particular study. We then perform a cut-and-count analysis 
on the event numbers in these enlarged signal regions. As input information we take into 
account the number of observed events, expected backgrounds, efficiencies, and systematic 
uncertainties, as published by ATLAS and CMS. Where necessary information is missing 
from the publications, we use estimations to the best of our knowledge. In App. [A] we list 
the data we use in more detail. 

For the actual fit, we follow a frequentist approach. In the absence of systematic uncer¬ 
tainties, the number of events observed in the signal region, n, follows a Poisson distribution 


p(n) 


H n exp (—p) 


( 1 ) 


Here p = b + a s ■ BR ■ L ■ e is the number of expected events, including the number of 
expected background events b, and the product of the signal production cross section a s , 
branching ratios BR, integrated luminosity L , and efficiency as well as acceptance factors e. 
The systematic uncertainties on the background prediction b and on the signal efficiency e 
are approximated by Gaussian distributions. These nuisance parameters are marginalized. 
We neglect correlations between the different systematic uncertainties. 

We then calculate a p value for each signal cross section in each different final state, 
using at least 10 000 pseudo-experiments. This gives us the significance of the deviation 
from the SM, best fit points, and confidence regions at 68 % and 90 % CL. We also calculate 
upper limits in a modified frequentist approach (CLs method |22j). In some of the signal 
regions, our simple cut-and-count strategy leads to model limits that are stronger than those 
presented by the experiments. In this case, the systematical background error and the signal 
efficiencies were rescaled to End agreement. Since many of our input parameters involve 
rough estimations, we have checked that the final results are stable under the variation of 
these input values. 

In a next step, we perform combined fits to all studies that contribute to a particular 
final state. For the combination of individual p values we use Fisher’s method. Again, we 
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present the results in terms of p 0 values, best fit points, and confidence regions at 68 % and 
90 % CL. 


II. 1 . Vector boson pair production 


WW resonance analyses 


Analysis 

Expected 
95% CLs [fb] 

Observed 
95% CLs [fb] 

Excess 

significance [a] 

Fitted cross 
section [fb] 

ATLAS hadronic [Tj 

11.1 

19.6 

1.2 

5.2 

CMS hadronic [H| 

12.2 

17.9 

1.0 

6.0 

ATLAS single lepton j6j 

6.4 

5.9 

0.0 

0.0 

CMS single lepton [5] 

7.2 

8.1 

0.3 

1.2 


Table I. Expected and observed limits, significance of any excess present, and the potential signal 
cross section for WW pair production at 1.8 TeV. Values may differ from those reported in the 
experimental papers since we integrate over signal regions in the vicinity of 1.8 TeV depending on 
the resolution and binning of the selected analysis. 


WZ resonance analyses 


Analysis 

Expected 
95% CLs [fb] 

Observed 
95% CLs [fb] 

Excess 

significance [cr] 

Fitted cross 
section [fb] 

ATLAS hadronic [I] 

14.2 

25.8 

1.3 

6.9 

CMS hadronic [9] 

11.9 

17.5 

1.0 

5.8 

ATLAS single lepton [6j 

13.2 

12.4 

0.0 

0.0 

CMS single lepton [5] 

14.9 

16.8 

0.3 

2.4 

ATLAS double lepton [7] 

13.8 

20.5 

0.3 

2.9 

CMS double lepton [5] 

14.4 

27.4 

1.5 

10.0 


Table II. Expected and observed limits, significance of any excess present, and the potential signal 
cross section for WZ production at 1.8 TeV. Values may differ from those reported in the exper¬ 
imental papers since we integrate over signal regions in the vicinity of 1.8 TeV depending on the 
resolution and binning of the selected analysis. 

The relevant searches targeting the resonant production of a pair of vector bosons ei¬ 
ther require both bosons to decay hadronically [I] [5] or one boson leptonically j5H7j. Un¬ 
fortunately, branching ratio suppression dictates that searches where both bosons decay 
leptonically are currently uncompetitive. 

In this set of analyses, the ATLAS fully hadronic search presently contains the largest 
single excess with 3.4 a reported by the experiment for an invariant mass of 2.0 TeV. Our cut- 
and-count analysis focuses on a slightly lower invariant mass window, leading to a reduced 
peak significance of 2.4 cr, see Tab. IH In order to fit the different final states we make use of 
the fact that the analysis has a mild discrimination between W and Z bosons based on the 
invariant mass of a reconstructed ‘fat jet’. This leads to a slight preference for either a WZ 
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ZZ resonance analyses 


Analysis 

Expected 
95% CLs [fb] 

Observed 
95% CLs [fb] 

Excess 

significance [a] 

Fitted cross 
section [fb] 

ATLAS hadronic [1] 

9.7 

25.2 

2.4 

8.1 

CMS hadronic [9] 

11.7 

17.1 

1.0 

5.7 

ATLAS double lepton [7j 

6.7 

10.0 

0.3 

1.4 

CMS double lepton [5] 

7.0 

13.4 

1.5 

4.9 


Table III. Expected and observed limits, significance of any excess present, and the potential signal 
cross section for Z pair production at 1.8 TeV. Values may differ from those reported in the 
experimental papers since we integrate over signal regions in the vicinity of 1.8 TeV depending on 
the resolution and binning of the selected analysis. 



Figure 1. p values for a 1.8 TeV resonance decaying into a combination of WZ and WW pairs 
(left), or ZZ and WW pairs (right). The 68% (90%) preferred region is shown with a solid (dotted) 
line. The star represents the best fit point. 


or ZZ final state to explain the excess present with a cross section of ~ 7-8 fb, see Tab. [Li 
and III Since a smaller peak is seen in the purely WW channel, a smaller cross section of 
~ 5 fb is found for this channel, see Tab. |TJ We again note that many of the events seen in 
the excess region are shared between all three final states; thus the total cross section in the 
excess region is substantially smaller than if we simply sum the three fitted cross sections 
together. 

The CMS fully hadronic analysis search is very similar and also finds an excess in the 
same mass range, although with a slightly smaller signal of ~ 1 a, see Tab. 0IH In this 
analysis no real preference is seen for any of the different final states and a cross section of 
~ 6 fb fits all three equally well. 

More discrimination of the final states is available by using the semi-leptonic searches. 
Here, the one-lepton analyses require that a W boson is present while the two-lepton searches 
reconstruct at least one Z boson. If we first examine the one-lepton searches we see that no 
excess is seen by ATLAS [6], while CMS [5] only has a very mild excess of ~ 0.3 a. These 
searches place a significant constraint on the WW final state with ATLAS and CMS giving 

















6 



o(pp -» X) x BR(X -> ZZ) [fb] o(pp -> X) x BR(X -> HW) [fb] 


Figure 2. p values for a 1.8 TeV resonance decaying into a combination of WZ and ZZ (left), or 
HZ and HW pairs (right). The 68% (90%) preferred region is shown with a solid (dotted) line. 
The star represents the best fit point. 


limits around 6-8 fb at 95% CLs. Removing one W for the WZ final state relaxes this limit 
to 12 fb, see Tab. |TT[ 

In the two-lepton case, CMS [5j sees an excess of 1.5 a which leads to a fitted cross section 
of 10 fb if we assume a pure WZ final state, see Tab. |H] or 5 fb for a pure ZZ final state, 
see Tab. |III| The ATLAS search has similar sensitivity [7] but only observes a very small 
excess of 0.3 a. 

In Fig. [I] and Fig. [2] (left panel) we show the combined cross-section fit to all channels. 
Fig. [ljshows that both WZ and ZZ final states fit the data well with a cross section of ~ 5 fb. 
However, we also see that a pure WW signal is disfavoured and could only describe the data 
in combination with another signal. The reason that the WW explanation is disfavoured is 
two-fold. First, the ATLAS and CMS single lepton analyses set an upper limit around 6 fb 
at 95% CLs, but a cross section of this magnitude is required to fit the hadronic excesses. In 
addition, the CMS dilepton search [5] has a small excess that this channel cannot explain. 

Finally the left panel of Fig. [2]shows that an equally good fit is possible with contributions 
from both a WZ and ZZ final state, i. e. requiring both a charged and neutral resonance 
with similar mass. In combination, the standard model is disfavoured by ~ 2 cr when all 
vector boson pair production channels are included. 


II.2. Associated vector-Higgs production 

The experimental searches for a resonance that decays into a HV final state are varied. 
ATLAS looks for a Higgs that produces a bb pair with the W or Z probed leptonically (££, 
u£ or vv) |8]. CMS has a similar search for H —>■ bb but only examines the leptonic W 
channel [0]. To probe vector bosons more generally, CMS has a fully hadronic search m 
but this has a limited ability to discriminate between W and Z. Finally there is a CMS 
search for H —>• t + t~, again with a hadronic reconstruction of the vector boson [TO] , 

Out of these searches, only the CMS study with a leptonic W displays an significant 
excess with ~ 2 a at 1.8-1.9 TeV. Interpreting this as a resonance that decays to HW leads 
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HW resonance analyses 


Analysis 

Expected Observed Excess 

95% CLs [fb] 95% CLs [fb] significance [a] 

Fitted cross 
section [fb] 

ATLAS bb + (ff, , vv) [8] 

33.0 

30.0 

0.0 

0.0 

CMS bb + v£ g] 

18.6 

44.4 

1.9 

15.8 

CMS t + t~ + hadronic vector [10] 

36.1 

36.1 

0.0 

0.0 

CMS hadronic Higgs [IT] 

12.5 

13.2 

0.1 

1.0 


Table IV. Expected and observed limits, significance of any excess present, and the potential signal 
cross section for associated production of a Higgs and a IV at 1.8 TeV. Values may differ from 
those reported in the experimental papers since we integrate over signal regions in the vicinity of 
1.8 TeV depending on the resolution and binning of the selected analysis. 

HZ resonance analyses 

Analysis Expected Observed Excess Fitted cross 

95% CLs [fb] 95% CLs [fb] Significance [er] section [fb] 


ATLAS bb + (M, ui , uu) [S] 

15.0 

14.0 

0.0 

0.0 

CMS t + t~ + hadronic vector [10] 

31.8 

31.8 

0.0 

0.1 

CMS hadronic Higgs [11] 

12.2 

12.9 

0.1 

1.0 


Table V. Expected and observed limits, significance of any excess present, and the potential signal 
cross section for associated production of a Higgs and a vector boson at 1.8 TeV. Values may differ 
from those reported in the experimental papers since we integrate over signal regions in the vicinity 
of 1.8 TeV depending on the resolution and binning of the selected analysis. 


to a fitted cross section of 16 fb, see Tab. IV However, the fully hadronic CMS search is 
slightly in tension with this result as it reports a limit of 13 fb on the same final state. In the 
analyses that are sensitive to HZ production, no significant excesses are seen, see Tab. |V} 
For this final state, the ATLAS semi-leptonic [8] and CMS fully hadronic m searches have 
similar sensitivities and set a 95% CLs limit of 14 fb and 13 fb respectively. 

Combining all these searches into a single fit, we plot the preferred cross sections for HW 
and HZ productions in the right panel of Fig. [2j We find that the best fit point has a cross 
section for HW production of ~ 5 fb, but is also compatible with the SM background at 
~ 1 a. Since the only excess is seen in a channel compatible with a W in the final state, 
there is no evidence in the data for a signal in the HZ channel. 


II.3. Dijet production 


Both the dijet search by ATLAS |3] and that by CMS [2j see an excess in the invariant 
mass distribution around 1.8 TeV. In our analysis, we find the excess in CMS is slightly 
more significant at 1.9 a compared to 1.5 a in ATLAS. However, since the CMS analysis is 
slightly more sensitive to a signal, the fitted cross section is actually smaller at ~ 90 fb, 
compared to ATLAS with ~100 fb, see Tab. VI In any case, the two signals seen by both 
experiments are remarkably similar and the combined best-fit signal cross section of ~ 90 fb 
can be seen in Fig. VlTJ 


















Dijet resonance analyses 


Analysis 

Expected 

Observed 

Excess 

Fitted cross 


95% CLs [fb] 

95% CLs [fb] 

Significance [c] 

section [fb] 

ATLAS dijet (3] 

131 

217 

1.5 

101 

CMS dijet [2] 

92 

173 

1.9 

90 


Table VI. Expected and observed limits, significance of any excess present, and the potential 
signal cross section for dijet production at 1.8 TeV. Values may differ from those reported in the 
experimental papers since we integrate over signal regions in the vicinity of 1.8 TeV depending on 
the resolution and binning of the selected analysis. 

II.4. Associated top-bottom production 


t — b resonance analyses 


Analysis 

Expected 
95% CLs [fb] 

Observed 
95% CLs [fb] 

Excess 

Significance [er] 

Fitted cross 
section [fb] 

ATLAS hadronic t [12] 

155 

203 

0.6 

31 

ATLAS leptonic t [13] 

138 

101 

0.0 

0 

CMS leptonic t [2] 

76 

67 

0.0 

0 


Table VII. Expected and observed limits, significance of any excess present, and the potential signal 
cross section for resonant production of a top and bottom quark at 1.8 TeV. Values may differ from 
those reported in the experimental papers since we integrate over signal regions in the vicinity of 
1.8 TeV depending on the resolution and binning of the selected analysis. 


The final analyses that we consider are the studies that look for the resonant production 
of a tb final state, see Tab. VII ATLAS has two searches that focus on this signature, one 
that looks for the hadronic decay of the t [12] and another that considers the leptonic channel 
[13] . On the CMS side, only one study exists, concentrating on the leptonic decay [13] , At 
the current time, only the ATLAS hadronic search contains a small excess in the region of 
1.8 TeV. However, since this search was expected to have the poorest sensitivity, it is likely 
that this is purely a statistical fluctuation. 

The strongest bound comes from the CMS leptonic search with an upper limit on the 
cross section of 70 fb at 95% CLs. However, we should also pay particular attention to the 
ATLAS leptonic search since this heavily influences our final model fits and interpretations. 
The reason is that at 1.8 TeV, the search records a signal 1.8 a smaller than expected. 
We hesitate to call this an ‘underfluctuation’ because the search systematically measures a 
cross section 2 a less than the background prediction across the whole mass probed (500- 
3000 GeV). Due to the construction of the CLs limit setting procedure, which weights the 
likelihood according to the agreement of the signal with the background, the systematically 
high background prediction does not result in a hugely significant shift in the 95% CLs limit 
(expected limit: 138 fb, observed limit: 101 fb). 

However, in our cross-section fit, the search has a far greater effect on the overall p value, 
since any additional signal predicted in this channel will be heavily penalized. For this 
reason, we perform two fits to the tb final state, one which includes the ATLAS leptonic 
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o(pp -> X) x BR(X -»jj) [fb] o(pp -> X) x BR(X -»jj) [fb] 

Figure 3. p values for a 1.8 TeV resonance decaying into a combination of jj and tb. The 68% (90%) 
preferred region is shown with a solid (dotted) line. The star represents the best fit point. Left: 
including the ATLAS leptonic tb search that displays a systematically ~1.8cr improved exclusion. 
Right: excluding the ATLAS leptonic tb search from the fit. 

search, shown in the left panel of Fig. [3| and one without, see the right panel of Fig. [3} 
In both plots we see that the best fit point is found without a signal present. However, 
including the leptonic ATLAS search results in a 1 a allowed cross section of ~ 30 fb but 
removing the search allows this to increase to ~ 60 fb in combination with the dijet result. 


II.5. Channel comparison 


Fitted cross sections 


Process 

Fitted cross 

Upper bound 


section [fb] 

(90% CL) 


pp —fX—t WZ 1 

r 7+3.6 

1 -3.3 

11.8 

pp — > X — > ZZ 1 

5.0±H 

11.3 

pp —> X —> WH 

4 5 + 5-2 
+•'-'- 4.0 

15.5 

pp -A X -A jj 

01+53 

45 

170 

pp —> X — > tb 

n+n 

u -0 

38 

pp —* X — > tb (without ATLAS bbiv [IB]) 

n+ 39 

u -0 

60 


Table VIII. Fitted cross sections to various final states found by combining all relevant channels. 
For the tb final state we present the results both with and without the ATLAS leptonic search m 
since this study contains a large under-fluctuation that significantly alters the result of our final 
fit. 1 The WZ and ZZ channels contain a significant overlap in signal regions and should not be 
considered as independent measurements. 

By comparing the various fitted cross sections we can try to provide guidance on the 
kind of model required to fit the current data. We also note that this is not completely 













10 


speculative since a combined analysis of the above searches finds that the standard model 
has a 2.9 a discrepancy with the data due to various excesses present. 

We first note that the required cross sections to correctly fit the data are very similar 
for WZ and ZZ (which can be considered as roughly the same measurement) and WH, see 
Tab. VIII In our opinion this seriously motivates a model in which the resonant production 


particle carries charge, since the same particle can then be responsible for both final states. 
In addition, the fact that awz ~ &wh suggests to us that a model that predicts an equal 
branching ratio to these two modes should be considered. 

In comparison, the dijet cross section for the signal is over an order of magnitude larger. 
In terms of finding a model to fit these excesses this is convenient, since the easiest way 
to produce such a high mass resonance will be through a quark coupling of appropriate 
strength. The same coupling will then automatically lead to a decay into a dijet final state. 

Combining the WZ, WH , and dijet final states thus naturally leads to a model with a 
charged resonance. Working with the principle of simplicity, we may expect the couplings 
of this resonance to be flavor diagonal and therefore predict that a t b ~ \&jj- Unfortunately 
there is no evidence for a signal in the tb final state and the 1 a preferred region only ext ends 
to ~ 10 fb, which is roughly However, we refer the reader to the discussion in Sec. II.4 

where we note that the ATLAS leptonic search finds a cross section that is systematically 
2 a below the background prediction. Consequently, the fit is heavily influenced by this 
single result and we believe that it is also wise to study the cross section measurement when 
this analysis is removed. In this case, the 1 a preferred region extends up to ~ 40 fb and is 
therefore perfectly compatible with the observed dijet signal. 


III. INTERPRETATION IN THE LEFT-RIGHT SYMMETRIC MODEL 


We now turn towards a specific model and interpret the observed excesses in the Left- 
Right Symmetric Model [19] ■ In Sec. III. 1 we summarize the structure and key phenomeno¬ 


logical properties of this framework. We then compare its predictions to the experimental 
data and perform a model parameter fit in Sec. III.2 Finally we discuss the prospect of 


observing signatures for this model during the upcoming LHC Run II in Sec. III.3 


III.l. Model essentials 

In the Left-Right Symmetric Model (LRM) framework [T!5], the SM electroweak gauge 
group is extended to SU(2) L x SU(2) R x U( 1)' with the usual left-handed (right-handed) 
fermions of the SM transforming as doublets under SU(2) L ( R y For instance, the left-handed 
(LH) quarks transform as (2,1,1/3), whereas the right-handed (RH) quarks transform as 
(1,2,1/3). The corresponding gauge couplings are denoted by ql,r , which are not in general 
equal, as well as g'. As a consequence of this L — R symmetry, a RH neutrino must necessarily 
be introduced for each generation; neutrinos are thus naturally massive in this framework. 

The SM is recovered by the breaking SU(2) R x U{1)' —> U(l)y, where additional Higgs 
fields are required to break this symmetry generating the masses of the new gauge holds W R , 
Z R . This breaking relates the g R and g' couplings to the usual SM hypercharge coupling 
r/y. Note that since gj, and gy are both well-measured quantities, only the ratio n = g R /gL 
remains as a free parameter in the model. The new Higgs fields responsible for SU(2) R 
breaking are usually assumed to transform either as a doublet, (1,2,—1), or as a triplet, 
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(1,3, —2), under the SU(2)r gauge group, and the neutral component obtains a multi-TeV 
scale vev, vr. The corresponding LH SU(2) L doublet or triplet Higgs fields, (2,1, — 1) or 
(3,1,—2), which must also be present to maintain the L — R symmetry, are assumed to 
obtain a vanishing or a tiny, phenomenologically irrelevant vev, which we will set to zero 
(i. e. vl = 0) below. 

The immediate impact of the choice of doublet versus triplet Higgs breaking is two-fold. 
In the triplet case, the RH neutrinos N e ^ r can obtain TeV-scale Majorana masses through 
the triplet vev, and thus see-saw-suppressed masses can be generated for the familiar LH 
neutrinos, which are now themselves necessarily of a Majorana character. Naively, we might 
expect the Wr and Ng masses to be of a similar magnitude. Furthermore, when triplet 
breaking is chosen we can easily identify U( 1)' with U(1)b~l- If, on the other hand, doublet 
breaking is assumed, then the LH and RH neutrinos must necessarily pair up to instead 
form Dirac fields. 

In addition to the many low-energy implications of this choice (such as potentially observ¬ 
able neutrinoless double beta decay in the Majorana neutrino case, or RH leptonic currents 
appearing in /j decay in the Dirac case), this selection has an immediate impact at colliders 
that is relevant for our analysis. If the neutrinos are Dirac, then the decay Wr —>■ 
occurs where the ug appear as missing energy. The LHC Run I searches for this mode con¬ 
strain the mass of Wr to lie beyond the range of interest for the present analysis, even 
for very small values of k ~ 0.15 (which he outside the LRM physical region as we will 
discuss below). This highly disfavors the doublet breaking scenario from our perspective. 
However, in the triplet breaking case, we instead find the leptonic decay now is of the form 
Wr —> £±Ng, so that the search reach depends on the relative ordering of the Wr and Ng 
masses. If the mass relation m N( > mw R is satisfied, then the Wr has no on-shell, two-body, 
leptonic decay modes and these LHC Run I search constraints are trivially avoided. If the 
RH neutrinos are below the Wr in mass, then ££jj final states will be produced. CMS [23] 
has observed a potential excess in this mode for the case £ = e but not for £ — ji. To 
interpret this as a real signal in the present scenario would require the various Ng to be 
non-degenerate, so that N e (IVjJ is lighter (heavier) than Wr [23TI28] . This, however, would 
naively lead to a predicted branching fraction for the process fi —>■ ey [29] which is far larger 
than the current experimental bound [30], unless the flavor and mass eigenstates of the Ng 
are extremely well-aligned. Furthermore, since the Ng are Majorana states, their decays 
produce the final states of like-sign as well as opposite sign leptons with equal probability, 
which is not what CMS apparently observes. To avoid these issues here we will thus assume 
both triplet SU(2)r breaking and that the relation rriRf e > rriw R is satisfied. 

Going further, we note that since the LH (RH) SM fermions transform as doublets under 
SU(2)i,(r), their masses must be generated by the introduction of one or more bi-doublet 
scalars, i. e. fields transforming as doublets under both SU( 2)l(r) groups simultaneously as 
(2,2,0). The vevs of these bi-doublets (with each bi-doublet, <Lj, having two distinct vevs, 
ki, 2 i which are of order of the electroweak scale) break the SM gauge group in the usual 
manner and act similarly to those that occur in Two-Higgs Doublet Models [3TJ. Given a 
sufficiently extensive set of bi-doublets, it is possible to construct models wherein the CKM 
matrices in the LH and RH sectors are uncorrelated. However, the relationship \V^\ = 
is the more conventional result if we want to avoid flavor-changing neutral currents; we 
will assume the validity of this relationship in the analysis below to greatly simplify the 
discussion as this additional parameter freedom is not needed here to explain the data. 
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We can now write the full W — Wr mass matrix in a generic manner as follows: 


M W *=( a m \ 

V Pwrriw m^ R ) 


( 2 ) 


Here one finds that = g\ Ei(^n + ^2j)/4, which we note is the would-be SM W mass, and 
correspondingly rriy VR = g R (2v R + + fcfj))/4. Note that, apart from a 0(1) coefficient 

, the off-diagonal term is proportional to rn 2 v . The reason for this is that the off-diagonal 
terms in the mass matrix are also generated by the vevs ki^i, so they are naturally of the 
order of the weak scale; one finds explicitly that 


2 E k H k 2l 

E Mi + k 2 2i )- 


(3) 


To diagonalize this matrix we rotate the original W, W R fields into the mass eigenstates W\p 
(where W\ is identihed as the well-known lighter state) via a mixing angle cj) w given by 


tan 2 <j) w 


-2j3 w m 2 w 


771 


W R 


m 


w 


(4) 


When rri 2 V/{ m as in the case under consideration, we obtain that (f> w ~ — f^w^w / m w R ) 2 ■ 
Note that in the most simple, single bi-doublet case, we fold that j3 w = 2Ktan/3/(l + tan 2 /3). 
Here we have defined the ratio of the kip vevs as tan/3, as usual. Although W \- 2 are the 
mass eigenstates, for clarity we will continue to refer to them as W, Wr. 

We can perform a similar analysis in the Z — Zr mixing case. This is simplified by first 
going to the basis where the massless photon is trivially decoupled, reducing the original 3x3 
mass matrix to one which is only 2x2. Then the SM Z couples as usual as —(T^l — %wQ), 
where c w = cos 9 W , Q is the electric charge, x w = s 2 ^ = sin 2 9 W and T^l is the 3 rd component 
of the LH weak isospin. Recalling Q = T^l + T^r + (B — L)/ 2 we can write the analogous 
Z R coupling as 


0 Zr = — [k 2 - (1 + k 2 )x u ] 1/2 [x w T 3L + k 2 (1 - x w )T 3R - x w Q}. 
Cw 


(5) 


Interestingly, since we know that v R A; 2 2i , the mass ratio of the physical Wr and Zr is 
given to a very good approximation by simply setting the k\ 2i —0, i. e. 


™ 2 z r = K 2 {l-x w )p R 
m 2 Wii k 2 (1-x w )-x w 


( 6 ) 


with the values of p R = 1(2) depending upon whether SU(2)r is broken by either Higgs 
doublets (or by triplets): in this work Pr = 2 follows from our assumption of triplet breaking. 
We demonstrate this relation between k and the physical masses in Fig. [4j 

In further analogy with the W — Wr case we fold that the Z — Z R mixing angle is given by 
c i> z — -ft(mz/m^) 2 when Zr is heavy. /3 Z is again an Oil) parameter which is generally a 
ratio of the various bi-doublet vevs. In the case of a single bi-doublet, the general expression 
for (3~ simplifies significantly to /3 Z = —[k 2 — (1 + k 2 )x w ] 1//2 . Note that for k values near 
the theoretical minimum (as we will discuss next) we fold that j3 z is very small, implying a 
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Figure 4. mz R as a function of k for m\y R = 1900 GeV. The physical region is shaded in grey. 


further suppression of Z — Zr mixing in this case. 

Examining the expressions for both the Zr/Wr mass ratio, as well as that for the Zr 
couplings, we see that k > [x w /(1 — Xw)] 1 ^ 2 ~ 0.55 is required for the fields to remain physical 
as alluded to in the discussion above. Below this value the Zr coupling becomes imaginary, 
see Eq. (J5]) , and is negative, c. f. Eq. (J6]) . This theoretical requirement will play an 
important role in the discussion of our fit results below. 

In this scenario, the partial width of the Wr into quark pairs is given by [32, 33] 


r(IE+ -> ud) = r(IE+ -> cs) = 3 k 2 A 
r(VE+ tb) = 3 k 2 A ^1 + 


x + ®s( m WR ) 

7r 



) 


2 m w R ) 


(7) 

( 8 ) 


where A = GF'niwmw R / (67ta/ 2) is an overall constant. 

Calculating the decays into diboson states is a little more complicated since correctly 
obtaining the effective WrWZ coupling in the LRM is subtle. As in the SM, the trilinear 
couplings of the gauge bosons arise from the non-abelian parts of the kinetic terms for the 
gauge fields, in particular, from the part of the covariant derivative containing the gauge 
fields acting on themselves. In the basis where the massless photon explicitly appears, the 
covariant derivative is given by 


V = d- ieQA - -^=g L T± ■ W T + (L —>■ R) — i^(T 3L - x w Q)Z - iO Zr Z r , (9) 

v 2 c w 

where we have suppressed the Lorentz index and where the Zr coupling operator Oz R is 
given in Eq. (|5j) . V acting on the W, Wr generates both a W + W~ Z coupling, as in the SM, 
as well as a W^W^Z coupling. Since W , Wr have Q — 1 and T 3L (W, Wr) = 1, 0, these two 
couplings are simply — (1 — ic^and ^r(—x w ), respectively. In terms of the mass eigenstates 
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W = W\c + IV 2 s and Wr = VV^c — W\s (where s = sin 4> Wl c = coscp w ), corresponding to 
the physical masses m lj2 , the off-diagonal WfW^Z coupling can be obtained by combining 
these two individual contributions. We obtain 

— [cs(l-x VJ ) + (-sc)(-x w )] = ( ^ L . (10) 

c w c w 

This reproduces the result obtained some years ago [3TJ. 

In the expressions above, we have ignored any Z — Zr mixing since it is numerically small 
in the parameter range of interest in our scenario since both \f3 z \ <C \/3 w \ and m^ VR /m 2 ZR 1. 
However, we note that by identifying Z with Z 1; an additional contribution will arise from 
the WfiWftZn coupling (the corresponding W + W~ Z R coupling is absent as can be seen 
from the structure of the operator Oz R ), but this interaction is relatively suppressed by an 
additional factor of the Z — Zr mixing angle, 0~, so we will ignore this term in our analysis. 

In the corresponding partial width for the decay Wr —» WZ (i. e. for W 2 —» W\Zi) 
the above coupling will appear quadratically and is always accompanied by an additional 
factor of m\/{m\m 2 z ) arising from the longitudinal parts of the corresponding gauge bo¬ 
son polarization appearing in the final state. Now, since (j) w ~ and mL = 

m w,w R [ 1 + 0( m w/ m w R )\i us i n g the SM relation mw = c w mz, we see that the large mass 
ratios will cancel, leaving us with just an overall dependence of 

In the single bi-doublet model this reduces further to ~ g 2 n sin 2 2/3 [1 + 0(m 2 v /rn 2 Vij )\ with 
tan (3 being the ratio of the two bi-doublet vevs as defined above. 

As shown in explicit detail recently in Ref. [35j (which we have verified), the corresponding 
square of the WrWH coupling, where H is to be identified with the (almost) SM Higgs, 
is given by ~ gj i cos 2 (a + (3), where a is the mixing angle of the two Higgs doublet model 
|31j . This is at the same level of approximation where higher order terms in the gauge boson 
mixings are neglected. Going to the Higgs alignment limit, i. e. a ~ (3 — tt/2, to recover the 
SM-like Higgs, one observes that cos 2 (a + (3) —> sin 2 2/3. This demonstrates the equality of 
the effective WrWH and WrWZ couplings up to higher order terms in the various mass 
ratios, as is required by the Goldstone boson equivalence theorem. 

To cut a long story short, the partial decay widths of the Wr into diboson states are 
given by 


r(w£ -> w + z) = 


A 2 


r, m 2 w + m 2 z + 


m 


w 


m%y 


m 


x 1 + 10- 


w R 

m'w + m 
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"Wr 
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( 11 ) 


( 12 ) 


with a w = (cs)(m^ /fl /m^ / ). an depends on the details of the Higgs sector; however, the 
equivalence theorem requires that clr = a w + 0( m w/ m w R )- 
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Figure 5. Fit of the LRM parameters k and sin (f)\v to the searches listed in Tab. VIII with 
rnwii = 1800 GeV. The region compatible with data at 68% (90%) CL is shown with a black solid 
(dotted) line. The white dashed line denotes the theoretical limit k > sw/cw — 0.55. The star 
represents the best fit point. Left: fit including the ATLAS leptonic tb search that systematically 
disfavors the full parameter space as well as the background-only hypothesis with ~1.8cx. Right: 
fit excluding the ATLAS leptonic tb search. The white solid line shows the 95% CLs limit from 
the ATLAS leptonic tb search (the region left of the line is allowed). 


III.2. Is that it? Fitting the Left-Right Symmetric Model to data 

We now investigate whether the excesses observed at the LHC can be explained by a 
Wn resonance of the Left-Right Symmetric Model. In section [TT] we concluded that the data 
favors a charged resonance with approximately equal branching ratios to WZ and WH. 
This agrees very well with the predictions of the LRM, where these branching ratios are 
equal up to corrections. The excess observed in the dijet channels, together 

with the not so restrictive bounds in the tb searches, are also promising. 

In such an interpretation, mw R should take on the value 1800 -1900 GeV to be compatible 
with the observed excesses. A mass of 2000 GeV, as favoured by the ATLAS diboson excess 
alone, appears to be disfavoured by constraints from the other channels, such as the semilep- 
tonic diboson searches. We will perform fits for the two scenarios mw R = 1800,1900 GeV, 
in the understanding that the limited number of events, width effects and experimental 
resolution makes it difficult to pin down the mass more precisely. 

In order to test the compatibility of the LRM with the data, we include the parameter 
space of the LRM in our cross-section fit. Again, we first calculate the compatibility of a 
parameter point with the observed event numbers in all individual analyses. This is followed 
by a combination of the results, giving an overall p value which we present the results in 
terms of best fit points and confidence regions at 68% and 90% CL in the parameter space. 

The narrow width approximation is used throughout. We fold that the width of the Wr is 
of order 1 - 2% of its mass in the best-fit region, making the error due to this approximation 
sub-dominant to the other uncertainties present. We calculate the production cross section 
of the Wr using the MMHT2014 NNLO pdfs [36] with constant NNLO k factors [37] while 
the branching ratios are calculated using Eq. 0-©. We do not assume an = a w , but 
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Figure 6. Fit of the LRM parameters k and sin bw to the searches listed in Tab. VIII with 
rnwii = 1900 GeV. The region compatible with data at 68% (90%) CL is shown with a black solid 
(dotted) line. The white dashed line denotes the theoretical limit k > sw/cw — 0.55. The star 
represents the best fit point. Left: fit including the ATLAS leptonic tb search that systematically 
disfavors the full parameter space as well as the background-only hypothesis with ~1.8cx. Right: 
fit excluding the ATLAS leptonic tb search. The 95% CLs limit from the ATLAS leptonic tb search 
is outside the plotted region and does not affect the preferred region. 


allow for corrections of order 0{m^ v /rri^ v ^) by introducing a parameter £ using 


(Xh 


m 


w 


771 


W R 


(13) 


£ is allowed to float in the range [—10,10] with a flat prior and is profiled over in our results. 

In Fig. [5] and [6] we present the results from our fit for m WR = 1800 GeV and mw R = 
1900 GeV respectively. In the left panels we show the overall agreement with data when all 
experimental studies presented in the previous section are included. In this case, no part 
of the parameter space of the LRM is compatible with data at 68% CL. This is due to the 
tension between the dijet excess in the ATLAS and CMS searches and the ATLAS tb search 
in the leptonic decay mode. However as we have argued in Sec. II.2 this ATLAS tb search 
finds a rate consistently 2 cr below the background expectation, and while this does not lead 
to a strong CLs limit, it severely punishes the overall p value in our fit for the full parameter 
space. 


We therefore also present results where this single ATLAS tb search is excluded from the 
fit itself, where we explicitly check that the best-fit regions are not excluded by the CLs 
limits from this study. The results are shown in the right panels of Fig. [5] and [6} There is 
now a well-defined region where the LRM agrees very well with all searches and can describe 
the observed excess while satisfying the constraints from the other searches. For mw R = 
1800 GeV this region is roughly given by 0.4 < k < 0.55 and 0.00011 < sin (f> w < 0.0015. 
For rriw R = 1900 GeV the smaller production cross section allows for larger couplings 
0.5 < k < 0.65 and 0.0012 < sin <%, < 0.0016. 

These preferred couplings fall into a special place in the parameter space: as described 
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Figure 7. Fit of the LRM parameters mz R and sin <f>w to the searches listed in Tab. VIII with 
rnwii = 1800 GeV (left) and mw R = 1900 GeV (right). The region compatible with data at 68% 
(90%) CL is shown with a black solid (dotted) line. The star represents the best fit point. We 
exclude the ATLAS leptonic tb search from the fit, but the white solid line shows the 95% CLs 
limit from this search (the regions right of the lines are allowed). 


in Sec. III. 1 the theory requires k > sw/cw — 0.55 to be consistent. This means that 
for mw R = 1800 GeV the preferred region at 1 a falls entirely in the unphysical regime. 
The situation is different in the rriw R = 1900 GeV scenario, where we find good agreement 
further from this boundary of the physically allowed region. 

Since the mass of the Zr is fixed by n and mw R ( as shown in Eq. (J6| and Fig. [4]) , we may 
also ask what mass the Zr must have for the LRM to be consistent with the observations. A 
coupling close to the boundary n > 0.55 corresponds to a very heavy Zr, while large values 


of k translate to lower Zr masses. In Fig. [7] we show the results of our fit, still excluding 
the ATLAS leptonic tb search, in terms of mz R and sin (j) w . 

Assuming mw R = 1900 GeV, the data permit a lower bound on mz R of around 4 TeV 
which is substantially above the masses probed so far at the LHC [38l 140] . mz R < 3 TeV. 
There is no upper bound on mz R from our fit. This directly follows from the fact that our 
preferred regions for k extend into the region k, ~ 0.55, where mz R becomes very large. 

Our fit also allows us to analyze the origin of the constraints. In Fig. [8] we show the 
constraints from diboson and fermionic final states separately. As expected from Eq. ([7]) 
and ([8]), the dijet and tb rates fix the overall coupling constant gR while the WZ and WH 
rates then set the mixing angle (j) w . We conclude this section with a comparison of the 
preferred regions for m WR = 1800 GeV and mw R = 1900 GeV in Fig. [9] 


III.3. Prospects for Run II 

As the LHC begins operations at 13 TeV, the prospects to discover or exclude the model 
presented here are excellent due to the steep rise in heavy particle production at higher 
energies. In this section we first estimate the amount of data that will be required to more 
thoroughly probe a possible Wr resonance at 1.9 TeV. We then more speculatively explain 
the prospects that the LHC may produce and detect the neutral Z R gauge boson. 
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Figure 8. Preferred regions at 68% CL for mw R = 1900 GeV. We separately show the constraints 
from bosonic (blue, dark hatched) and fermionic (green, light hatched) final states. The ATLAS 
leptonic tb search is excluded from the fit. Left: parameterisation with the LRM parameters n 
and sin^w The grey shaded region marks the theoretically allowed region of k > sw/cw — 0.55. 
Right: parameterisation in terms of mz R and sin^y. 



Figure 9. Preferred regions at 68% CL for two different values for the Wr mass. The ATLAS 
leptonic tb search is excluded from the fit. The grey shaded region marks the theoretically allowed 
region of k > sw/cw — 0.55. 
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At the 13 TeV LHC, the production cross-section for a Wr resonance at 1.9 TeV is over 6 
times higher than at 8 TeV. Consequently, if we assume that the background scales roughly 
in proportion with the signal, a mere 5 fb -1 will already probe the model in more detail 
than the current data set. Indeed, if no signal is observed, the dijet resonance search can 
already be expected to exclude our best-fit point at 95% CLs. This would also place the 
whole model under significant strain since it is the jj cross-section measurement that drives 
the overall coupling determination in our fit. If we do not see a continued excess here, the 
model is driven to couplings of unphysically small sizes for this value of the Wr mass. 

As Run II accumulates 10 fb -1 , a signal should be observed in the tb final state otherwise 
the model assumption of flavour diagonal couplings will start to be under significant tension. 
On the gauge coupling side, our best fit model point can easily be excluded by both WZ 
and WH searches with less than 15 fb _1 . 

A categorical 5cr discovery does require larger data sets. Perhaps surprisingly given that 
most of the theoretical excitement has revolved around the diboson excesses, in our model 
we can expect the jj final state to be discovered first. Indeed, using the current search as a 
baseline, we expect a 5cr discovery to be made with approximately 20 fb -1 if the current best- 
fit point is close to reality. A discovery in the tb final state would follow shortly afterwards 
with 30 fb -1 . Again the gauge boson final states require more data with approximately 
50 fb -1 required for confirmation of the WZ final state, while over 100 fb -1 is expected to 
be required before WH is definitively seen. 

More speculative is the discovery potential of the Zr resonance since our model best fit 
point lies so close to a critical theory region. As explained in Sec. III.1[ model consistency 
requires that k > 0.55, and depending on the mass assumed for the Wr resonance, the 
best fit n may be below this, see Fig. [5] and [6j The problem is that at as we head towards 
k —> 0.55, the mass of the Zr rapidly increases to a value far above the LHC collision energy. 

Nevertheless, there is still a region of the 1 a preferred fit value that allows for LHC 
relevant Zr masses, especially if we assume mw R = 1900 GeV, see Fig. [9] A scenario with 
mz R > 4000 GeV is entirely possible and in Fig. [lOjwe give an estimation of the Zr —* l + l~ 
cross section at yfs = 13 TeV as a function of mz R - We find that the relevant region begins 
already to be probed once the LHC collects 20 fb -1 at 13 TeV. With 100 fb -1 , the ~ 95% CLs 
exclusion region stretches to ~ 4500 GeV. 


IV. A CONNECTION TO DARK MATTER? 

An immediate question that arises when hints for new physics are found at the LHC is 
whether these hints could be connected to the physics of dark matter. In the context of the 
anomalies discussed in this paper, such a connection is not obvious. We will discuss four 
different scenarios in the following. 

a ) W.R-mediated DM interactions with SM partners. If thermally produced DM 
particles x° are coupled to the SM sector through a Wr resonance, they need a charged 
partner which they can annihilate into. x + could be a SM lepton, in particular the r, 
but in this case the DM mass would have to be smaller than m T to forbid a fast DM decay 
via x° ~ t t + (Wr —> qq) (the coupling of the Wr to quarks is needed to explain the LHC 
anomalies discussed in the previous sections). In the context of the Left-Right Symmetric 
Model, x° could for instance be identified with the third generation right-handed neutrino. 
Then, however, mixing between y° and the other right-handed neutrinos must be forbidden 
or strongly suppressed to avoid fast decays into electrons and muons. Moreover, the Wr- 
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M Zr (TeV) 


Figure 10. Zr production cross section times branching ratio into the e + e~ and p + pX states 
combined for y/s = 13 TeV, assuming mw R = 1900 GeV and 0.56 < k < 0.8, see Eq. Q. The 
dashed lines show the cross sections that predict 3 expected signal events after 20 (blue) and 100 
(green) fb _1 of data taking. Due to the low backgrounds, these lines indicate ~ 95% CLs exclusions 
if no signal is seen. 


mediated annihilation cross section, which is of order 

2 

777/ 

(crv re i) (x°X° t + t~) ~ —£ - 1(T 30 cm 3 /sec, (14) 

m w r 

would be far below the generic value for a GeV thermal relic, (<7f re i) ~ 5 x 1CU 26 cm 3 /sec BH- 
In order to obtain the correct DM abundance in this case, significant entropy would need 
to be produced in the early Universe after DM freeze-out. While this is certainly possible, 
it would require a dark sector with a larger particle content and much richer dynamics 
than envisioned here. Alternatively, RH neutrino DM could be produced via a freeze-in 
mechanism. A scenario along these lines has been studied in Ref. [42] • 

b) WR-mediated DM interactions with charged partners beyond the SM. If both 
the x° an d X + are new particles, a thermal freeze-out scenario is possible if the y + is < 10% 
heavier than the y°. Under this assumption, DM decay is forbidden, but freeze-out through 
X°~X + coannihilation in the early Universe is possible. Regarding the DM phenomenology 
today, neither direct detection nor indirect searches are expected to yield signals in this 
coannihilation scenario. At the LHC, however, DM could be detected in anomalous decays 
of the Wr if m x o + m x + < m\y R ■ The process of interest is pp —> Wr —> x° + (x ± —• y X°QQ')j 
where the 3-body x ± decay proceeds through an off-shell Wr. Possible search channels for 
this process are thus jets and missing transverse energy (MET), or the associated production 
of a top and bottom quark with MET. Both signatures are plagued by large SM backgrounds, 
mainly from vector boson + jets production, and it is hard to estimate their discovery 
potential without running full simulations. In addition, x + X~ P a i r production would lead to 
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final states involving multiple jets and MET, similar to typical signatures of supersymmetric 
models. Note that in this scenario the branching ratios of the Wr into SM particles would be 
reduced compared to our assumptions in the previous sections, potentially allowing slightly 
larger values of k to be consistent with data. 


c) Z- and Zr -mediated DM interactions. Since W' bosons typically come with 
neutral Z'-like partners, it is also interesting to consider Z'-mediated DM-SM interactions. 
We will do this in the context of the LRM discussed in Sec. Ill, but our conclusions easily 
generalize to other models, in particular to scenarios in which the LHC diboson anomaly is 
interpreted as being directly due to a Z' resonance. The DM candidate could again be one 
of the standard right-handed neutrinos Nr. In this case, however, we would again face the 
problem of fast DM decay through the Wr. Therefore, let us consider a scenario where a 
new fermion multiplet x with quantum number (1, 2, —1) under SU( 2) l x SU (2 )r x U (1) b-l 
is added to the model. The upper (neutral) component of this doublet is the DM candidate 
X°, the lower (charged) component y - is assumed to be sufficiently heavy for coannihilations 
to be negligible. Since the multiplet (x°,X~) has the same quantum numbers as the right- 
handed leptons, there is the possibility of an undesirable mixing between x° and Nr, which 
would reintroduce IT^-mediated DM decay. Such mixing can be forbidden by introducing an 
additional symmetry, for instance a dark sector Z 2 symmetry. Note that such a symmetry 
does not forbid a Majorana mass term for the y 0 , unlike for instance a U( 1) symmetry. This 
is crucial, because if y° was a Dirac fermion, its vector couplings to the Zr and (through 
gauge boson mixing) the Z would bring it into blatant conflict with direct detection results. 
For a Majorana DM particle, however, the dominant effect in direct detection experiments 
is spin-dependent scattering through axial vector interactions, for which limits are much 
weaker. 


Annihilation of the y° in the early Universe proceeds through s-channel Z and Zr ex¬ 
change into fermionic final states and into W + W~ , with the fermionic final state dominating. 
This implies in particular that for m x ~ mzj 2 or m x ~ mz R /2 the total annihilation cross 
section (cru re i) tot is resonantly enhanced. Close to the Zr resonance, the annihilation cross 
section to massless fermions is given by 
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( 15 ) 


Here g'^ v , g'j A and g' xA are the vector and axial vector couplings of the final state fermions 
and the axial vector coupling of the DM particle to the Zr, respectively. Explicit expressions 
for them have been given in Sec. Ill [32, [33]. Note that x° does not have vector couplings 


because it is a Majorana fermion. This the reason for the v^ el suppression in Eq. (15). It can 


be understood by noting that due to the Pauli exclusion principle, the two incoming DM 
particles can only be in an s -wave state if their spins are opposite. The final state fermions 
are, however, produced in a spin-1 state due to the chirality structure of the gauge boson 
couplings. Thus, either one of them has to experience a helicity flip (which is only possible 
for m,f ^ 0), or the initial state DM particles have to be in a p-wave state. Note that in 
the resonance region, the p-wave contribution proportional to V rel is dominant. The reason 
is that, on resonance, an on-shell, spin-1 Zr boson is produced, and this requires the DM 
particles to be in a spin-1 state as well. Outside the resonance region, the p-wave terms 
dominate in the early Universe, where (^rel) rv_/ 0.24, while today, where ( V rel) 10 6 in the 
Milky Way, it is the helicity-suppressed terms that give the main contribution. 
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This implies that the annihilation cross section today is several orders of magnitude below 
the thermal relic value, making indirect DM detection in this scenario extremely challenging. 

To compute (crn re i) tot , we have used FeynCalc [33] to evaluate the annihilation cross 
sections for the // and W + W~ final states. Note that, in doing so, we need not only 


the coupling constants appearing in the simplified expression Eq. (15), but also the DM 


coupling to the SM-like Z boson. It is given by its coupling to the Z R , multiplied by the 
Z-Z R mixing angle /3 z (mz/mz R ) 2 as discussed in Sec. III. The ZW + W~ coupling is at its 
SM value, while the Z R W + W~ coupling is again suppressed by a factor /3 z (mz/mz R ) 2 ■ 
Note that in evaluating the cross section for XX ~t W + W~, we include only the transverse 
polarization states of the internal and external gauge bosons. Thus we avoid having to 
include model-dependent diagrams with Higgs boson exchange. Since annihilation to W + W~ 
is subdominant by a large margin compared to annihilation to fermions, this approximation 
will not affect our results. 

We plot (cr^ re i) tot as a function of the DM mass in Fig. 11 and compare it with the 
value required for a thermal relic [4TJ. Here we assume the conditions at DM freeze-out; in 


particular, we take (rr el ) ~ 0.24 for the average relative velocity of the two annihilating DM 
particles. We End that even at the Z resonance, the helicity and velocity suppression leads 
to annihilation cross section several orders of magnitude below the required value for the 
correct relic density. Only if the DM mass is close to mz R l 2 mass, the resonant enhancement 
is large enough to make the annihilation cross section compatible with the observed relic 
density. On the other hand, this implies that, if the Zr boson in this scenario is indeed 
responsible for the coupling of DM to the SM, the model provides a strong indication for the 
value of m x . Let us remark again that also an annihilation cross section somewhat below 
the naive thermal relic value may be acceptable if the LIniverse goes through a phase of 
extra entropy production after DM freeze-out, thus diluting the DM density. 

It is also important to consider the dark matter-nucleon scattering cross section probed 
by direct detection experiments. Since y° is a Majorana particle with only axial vector 
couplings, the scattering will be spin-dependent. The cross section is 


3 m 2 ™ 2 


a X N - 


N m x 


7r(rriN + m x ) 2 


E 


L q=u,d,s 


A (9qAg X A , 9qA9 x A\ 


1 2 


(16) 


where N = p or n for scattering on protons and neutrons, respectively. As before, primed 
coupling constants denote couplings to the Z R and unprimed ones indicate couplings to the 
Z. The hadronic form factors A N q are taken from Ref. [44]. In Fig. 12 we show the cross 
section of DM interactions with protons as a function of the DM mass and of k and compare 
it to limits from the XENON-100 [45] and PICO-2L [46] experiments. Independent of the 
DM mass, we hnd that our scenario predicts a cross section a few orders of magnitude below 
current exclusion bounds. 


We can conclude that this simple scenario gives the correct relic density for DM close to 
the Z R mass, while being fully consistent with current limits from direct detection experi¬ 
ments. Unfortunately, it is also not in reach of these experiments in the foreseeable future. 
We also expect the y° to be too heavy to be within the discovery reach of the LHC. 

d) Minimal Left-Right Dark Matter. Another possibility for DM in the LRM is to 
introduce a pair of chiral fermion multiplets 


Xl~(3, 1 , 0 ), Xr~ ( 1 . 3 , 0 ), 


( 17 ) 
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Figure 11. The total DM annihilation cross section (<7u re i) tot as a function of the DM mass for 
different values of k = qr/ql- We include fermionic final states as well as annihilation to W + W~. 
Due to the velocity and helicity-suppression of the annihilation cross section, we find that only for 
m x ~ m2 JJ /2, the cross section can be large enough to avoid DM overproduction in the simplest 
thermal freeze-out scenarios. 
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Figure 12. The spin-dependent DM-proton scattering cross section as a function of the DM mass 
for different values of n = gn/gLi compared to spin-dependent direct detection limits from the 
XENON-100 [35] and PICO-2L [4S| experiments. 
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which share a common Majorana mass M due to left-right exchange symmetry and whose 
neutral member (s) can be identified as DM. This scenario was recently considered in Ref. [47] . 
As was discussed there, such a scenario actually leads to a two-component picture for DM. 
Prior to clectroweak radiative corrections, all members of the x multiplets are degenerate. 
An issue that arises for the right-handed component is the sign of these corrections, i.e., 
whether or not they drive the masses of the charged states below that of that of the neutral 
member for some choice of the parameter ranges, in particular, the value of M. Ref. m 
showed that for M\y lt = 2 TeV and assuming k = 1 one must have M below ~ 1.8 TeV, 
otherwise this mass splitting goes negative thus preventing us from identifying the neutral 
component as (part of) DM. It thus behooves us to determine if this result is robust when 
we lower the value of k to our range of interest. Employing the Wr and Zr couplings given 
above (as well as the general Zr — Wr mass relationship), we find that the relevant mass 
splitting is now given by the expression 


AM = q2 GrM 2 , 


M 


2y/2- 


7T 


K 2 f( r Wr 


K, 2 — (1 + K 2 ) 




f( r Z E 


x: 


Xn 


-f(r z ) ~ x w f{r 0 


where r* = Mi/M and f(r) is given by the integral 


. (18) 


dx 2(1 + x ) log[x 2 + (1 — x)r 2 ], 


(19) 


which we can evaluate analytically^ Fig. 13 shows the values we obtain for AM vs. M as 
we vary k over the relevant range as well as for the case of k — 1 for comparison purposes 
always taking Mw R = 1.9 TeV. Here we see that the mass splitting is always positive for the 
range of n values of interest. We also observe that the magnitude of AM decreases as k value 
increases. We note that at large masses, as k increases, the curves begin to bend downward 
with an ever increasing slope. For the mass range shown here AM always remains positive 
until values of ft ~ 0.88 are reached. 

e) DM in supersymmetric Left-Right models. A fifth possibility of introducing 
dark matter in the framework discussed here is to consider supersymmetric grand unified 
models based on Left-Right symmetry. In this case, the lightest neutralino is an excellent 
DM candidate, as has been extensively studied [48j |49j. However, its phenomenology is very 
similar to that of the lightest neutralino in the minimal supersymmetric standard model 
(MSSM), with no direct connection to the excesses seen at ATLAS and CMS. 


V. CONCLUSIONS 

In this study we have analyzed the observed resonant excesses that appear in different 
ATLAS and CMS search channels in the invariant mass region of 1.8-2.0 TeV. The most 
prominent of these displays a 3.4cr excess in the search for the hadronic decay of a WZ final 
state, though this peak is also sensitive to other diboson states. Both ATLAS and CMS 
End excesses in the dijet distributions at approximately the same mass, with a significance 
around 2.2a and la. There are further potential signals in CMS searches for semileptonic 
decays of vector boson pairs as well as for resonances decaying to a HW pair, significant 


1 It is interesting to observe that the mass splitting vanishes in the case when all the mass functions, /, are 
equal as it does in the case of the SM. 
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M (GeV) 

Figure 13. For the minimal left-right DM scenario, we show the radiative mass splitting AM 
between the charged and neutral components of the new nmltiplets xl and \R as a function of the 
tree level Majorana mass M. The different colored curves correspond, from top to bottom, to k 
values going from 0.56 to 0.71 in steps of 0.03. We always assume M\y R = 1.9 TeV. The result for 
ac = 1 is shown as the dotted curve for comparison. 


at 1.5 cr and 2.2 a. We have investigated potential scenarios that can explain these excesses 
while being consistent with all constraints from other searches. 

In the first part of onr analysis, we have performed a model-independent fit of the cross 
sections corresponding to the observed and expected event numbers in all searches from 
ATLAS and CMS that are sensitive to these and closely related final states, including the 
analyses that agree with background expectations. Our fit finds an overall tension between 
the SM and the data equivalent to 2.9 a. To explain the excess observed in vector boson 
pair production and associated Higgs production, a charged resonance is favoured over a 
neutral resonance. In particular, states decaying into W + W~ pairs are strongly constrained 
from semileptonic searches. The best agreement is found for a charged resonance with ap¬ 
proximately equal branching ratios into WZ and WH pairs, with fitted signal cross sections 
around 5 fb in each final state. The excess in the dijet distributions suggest a branching 
ratio into quarks or gluons that is larger by a factor of at least 10, which is welcome news, as 
a sizable coupling of the heavy resonance to quarks or gluons is necessary for a sufficiently 
large production cross section at the LHC. Resonance searches in the tb state do not observe 
an excess, but the upper limit still allows a decay into this channel of approximately half 
the size of the dijet branching fraction. 

As a next step, we have interpreted these signatures in the context of the Left-Right 
Symmetric Model based on the extended gauge group SU(2) L x SU(2) R x U{1)' as the 
resonance production of a new heavy charged gauge boson, Wr. Fitting this model to the 
data, we have found that a Wr of 1900 GeV is in good agreement with all analyses, if 
the right-handed coupling is in the range 0.35 < g R < 0.45 and the mixing between the 
W and Wr is of the order siiuj) w ~ 0.0015. This preferred region can be translated into 
mass constraints for the associated heavy neutral gauge boson Zr , for which we find a lower 
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bound of approximately mz R > 4 TeV. For a lighter W R , the bounds from the fit become 
stronger, requiring a smaller coupling g R and a heavier Z R . 

In the upcoming 13 TeV run, the LHC will be able to probe this potential signal with very 
little data. Already with an integrated luminosity of 5 fb _1 , the experiments should be able 
to exclude the dijet signal of our best-fit scenario, followed by sensitivity in the tb channel 
shortly thereafter and then in the diboson channels with statistics of roughly 15 fb _1 . For 
a 5cr discovery, we estimate that a luminosity of 20 fb _1 is needed in the jj channel, while 
the tb and especially the WZ and WH states require more statistics. Whether the 13 TeV 
LHC will be able to produce the Z R crucially depends on the value of n (or equivalently 
mz R )- The LHC will begin to probe the interesting parameter space for the new neutral 
gauge boson with integrated luminosities around 20 fb 1 , but it is possible that this gauge 
boson is too heavy to be accessible at the 13 TeV LHC. 

In addition, we have analyzed if this model can simultaneously explain dark matter. If 
the DM annihilates primarily through the W R , we require two new states, one charged and 
one neutral (the DM candidate). These would have to be relatively close in mass in order to 
be able to coannihilate and produce the correct relic density. In the case of a Z R mediator, 
we introduce a new fermionic doublet with a neutral component to act as the dark matter. 
Here we find that a resonant mechanism is required to produce the correct amount of dark 
matter and thus if this solution is realized in nature, we predict the mass ~ \mz R - 

The fact that a number of excesses in different search channels across two experiments 
can be explained by an existing, simple (and, some might argue, natural) extension of the 
Standard Model is exciting. After the first months of data taking at the LHC at 13 TeV we 
will know more. We eagerly await the discovery of symmetry restoration in the upcoming 
operations of the 13 TeV LHC! 


NOTE 

While this study was in the last stages of preparation, several other studies appeared 
that analyzed similar models [271 [501153] . In addition, other models have been put forward 
to explain the various excesses mm- 
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Fit input data 


Analysis 

Selection 

Mass bins [GeV] 

Obs. 

Bkg. 

(unc.) 

Eff. 

(unc.) 

ATLAS VV hadronic 0 

WW selection 

1750-2050 

13 

8.5 

1.3 

0.10 

0.04 

ATLAS VV hadronic [l] 

ZZ selection 

1750-2050 

9 

3.0 

0.8 

0.08 

0.02 

ATLAS VV hadronic [[] 

WZ selection 

1750-2050 

18 

10.0 

1.5 

0.09 

0.03 

CMS VV hadronic [9] 

Double tagged 

1780-2030 

108 

96.4 

5.0 

0.22 

0.04 

ATLAS VV, single lepton [6] 

Merged region 

1700-2000 

8 

9.1 

5.2 

0.27 

0.01 

CMS VV, single lepton [5] 

High purity 

1700-2000 

12 

12.3 

5.3 

0.26 

0.03 

ATLAS VV, double lepton [7] 

Merged region 

1680-2060 

1 

0.5 

0.1 

0.24 

0.03 

CMS VV, double lepton [5] 

High purity 

1700-2000 

7 

3.5 

0.4 

0.41 

0.06 

CMS VH bb + vt g] 


1700-2000 

3 

0.5 

0.4 

0.06 

0.01 

CMS VH —» t+t~ + hadronic vector lit 


1500-2000 

8 

8.3 

3.5 

0.37 

0.05 

CMS VH, hadronic Higgs [llj 

bb selection 

1690-2030 

28 

27.1 

4.1 

0.16 

0.03 

ATLAS dijet [3] 


1706-2030 

38326 

37998 

90.0 

0.16 

0.02 

CMS dijet (2) 


1678-1945 

114117 

113438 

100.0 

0.38 

0.04 

ATLAS tb, hadronic t [12) 

Double tagged 

1600-2000 

432 

410.6 

28.0 

0.05 

0.02 

ATLAS tb, leptonic t [13] 


1600-2000 

14 

31.5 

16.9 

0.06 

0.02 

CMS tb, leptonic t [14] 


1500-2000 

178 

187 

20.0 

0.13 

0.01 


Table IX. Input data into our fit. We give the observed and expected background events, the 
uncertainty on the background events, the product of all relevant acceptance and efficiency factors, 
as well as the uncertainty on this number. 


Appendix A: Fit input data 


In Sec. [TT] we have described a cross-section fit to all ATLAS and CMS analyses sensitive 
to diboson, VH, dijet and tb final states. This fit is based on a cut-and-count analysis in 
the mass bins around 1800 GeV. In Tbl. |IX| we give the input information our fit uses in 


more detail, including the exact selections and mass bins that were used and the number of 
observed and expected background events in these analyses. Where available, these numbers 
are based on the papers, conference notes, supplementary material, and HepData entries 
published by ATLAS and CMS. Missing pieces of information were roughly estimated based 
on the available information. In some cases the estimated signal efficiencies and uncertainties 
were rescaled to give better agreement with the limits published by ATLAS and CMS, but 
we checked that this modification does not affect the overall fit results significiantly. 


In analyses involving hadronic decays of gauge bosons, in particular in the ATLAS diboson 
search p], IF bosons can be reconstructed in the Z selection and vice versa. Based on 
Fig. 1 c) in [Ij we estimate this spill factor to be 


p(W reconstructed as Z) ^ p(Z reconstructed as W) ^ 
p(W reconstructed as IF) p(Z reconstructed as Z ) 


(Al) 


The WW, WZ and ZZ selections in [T] are not orthogonal. In fact, in the enlarged signal 
region there are only 2 observed events that are tagged in either of the WW or ZZ selections, 
but not in the WZ category (see Fig. 13 of the auxiliary material published with pQ). This 
is why for the statistical combination of different searches we follow a conservative approach 
and only include the results from the WZ selection. 
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